Thermomagnonic spin transfer and Peltier effects in insulating 
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Abstract - We study the coupled magnon energy transport and collective magnetization dynam- 
ics in ferromagnets with magnetic textures. By constructing a phenomenological theory based on 
irreversible thermodynamics, we describe motion of domain walls by thermal gradients and gen- 
eration of heat flows by magnetization dynamics. From microscopic description based on magnon 
kinetics, we estimate the transport coefficients and analyze the feasibility of energy-related appli- 
cations in insulating ferromagnets, such as yttrium iron garnet and europium oxide. 
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Introduction. — Most electronic devices rely on 
charge current flows controlled by applied voltages. It has 
been realized that spin-polarized charge flows can also be 
induced by fictitious electromagnetic fields due to mag- 
netic texture dynamics [lp]. Recently, the possibility of 
employing spin in electronic logic devices has also been 
suggested |3|. However, since it may be energy-costly 
to create and maintain spin imbalances electrically, the 
fictitious electromagnetic fields induced by magnetic tex- 
tures may become useful for such applications. Reciprocal 
motion of magnetic textures due to weak charge-current- 
induced spin transfer is governed by dissipative torques 
(so-called "/? terms") [4[j8] associated with microscopic 
spin misalignments, which will also play an important role 
in our study. 

The spin flows are known to coexist with heat flows 
in transition metals, where both fluxes can be effectively 
driven by the moving magnetic texture (9l[l0] . Similar 
effects have also been studied in magnetic semiconduc- 
tors [IT]. The reciprocal action of the spin-transfer torque 
(STT) on the magnetization [12] is extremely important 
for applications, e.g., spin-transfer torque memory and 
nonvolatile logic. Spin torques due to pure spin currents 
have been observed in experiments on spin pumping by 
magnetic precessi on [l3] or spin Hall effect along magnetic 
insulator surface |[14|. Thermal spin torques have been 
studied in magnetic nanopillars [l5] with recent develop- 
ments suggesting that they can be more efficient compared 
to electrically generated spin torques |16 . 



Spin and heat currents can be also induced in insulat- 
ing systems, such as yttrium iron garnet (YIG), by the 
external microwave magnetic field |14|. Alternatively, the 
temperature gradient can lead to a spin imbalance by the 
spin Seebeck effect |17J. The heat currents accompanying 
the spin-wave (magnon) flows in an insulating ferromag- 
net LU2V2O7 have been used to detect the magnon Hall 
effect |18|. Here, we are interested in the possibility to 
use fictitious electromagnetic fields induced by magnetic 
textures to control magnon spin flows in a manner similar 
to spin-polarized charge currents. 

In this Letter, we study interplay between magnon- 
carried heat and spin currents and magnetic texture dy- 
namics in ferromagnets by formulating a hydrodynamic 
phenomenological description of magnonic and thermal 
currents coupled to geometric gauge fields (the charge cur- 
rents are ignored in our description but can be readily 
reintroduced in case of metals |9|). By considering dissi- 
pative corrections ("/? terms"), we show that they play an 
important role as they do in electronic systems, enabling 
domain- wall (DW) motion below the Walker breakdown 
and heat pumping |9|. To justify our phenomenology, we 
also formulate a ground-up kinetic description of thermal 
magnons, which allows us to identify the phenomenolog- 
ical coefficients entering the hydrodynamic theory. The 
DW motion and heat pumping are studied in YIG and 
EuO and the feasibility of energy-related applications is 
discussed. 



p-1 



Alexey A. Kovalev 1 Yaroslav Tserkovnyak 2 




Fig. 1: The magnon current induced by temperature gradient 
exerts spin torque on the magnetization, according to the con- 
servation of angular momentum, which can lead to a domain- 
wall motion. The inverse effect of magnon current induced by 
magnetization motion is also possible, wherein the collective 
magnetic texture m(x,t) is controlled by effective field H. 



Phenomenological description of magnonic and 
thermal currents in textured magnets. — We con- 
struct a general phenomenological description of magnonic 
and thermal currents in textured magnets (Fig. [T]) based 
on several thermodynamic variables such as the direction 
of the slow (averaged over the magnonic excitations) spin 
density m s (the label s is dropped in this section), den- 
sity of magnons p and density of energy pu . The ferro- 
magnet can be taken out of equilibrium by applying tem- 
perature, chemical potential and magnetization gradients, 
while the equilibrium state can be topologically nontriv- 
ial, e.g., a magnetic DW or vortex. In the absence of bias, 
the ferromagnet would then evolve back towards equilib- 
rium according to the equations of motion. Following the 
standard route, we identify thermodynamic variables x k 
and their conjugates (generalized forces) X k = dS>/dx k by 
writing the entropy 8 production in the following form: 



(i) 



The local conservation laws of energy and (approximately) 
magnon number provide us with continuity relations: p = 
— Vj — (p — po)/r a and pu = — Vjjy, where we introduced 
the magnon- number (j) and energy (j^y) current densities, 
and r a corresponds to the life time of magnons (that are 
in excess of a local-equilibrium value po). We now write 
the rate of the entropy production in a standard manner 
1 191: 



S = 



/ 
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Vjc/ + lip + U ■ m 
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where the conjugate force corresponding to the magnetic 
(spin-density) direction m is denned as — £m§|j[/(j)=o = 
H/T. By straightforward manipulation, in which we in- 
troduce the modified energy current } q — }u — /ij, we arrive 
at the following equation for the rate of the entropy pro- 
duction: 



/ 



_ . n . v/i. n : 

V ( - I j q - — j - — • m 



(3) 



where in Eq. Q we integrated the term involving } q 
by parts, used the continuity equations and disregarded 



magnon decay oc (p — po)/r a . The latter is justified when 
either (i) the number of magnons is (to a good approxi- 
mation) conserved, which means that the relevant size of 
the system is smaller than the magnon decay length cor- 
responding to Tq,, or (ii) no build-up of magnons takes 
place due to fast relaxation or uniform current genera- 
tion (thus p ~ po). The forces conjugate to the fluxes can 
be immediately identified as — <9j g §| m j=o = — V (1/T) and 
— <9jS| m j g =o = V/i/T. Formally, Eq. ^ is identical to the 
one used in Ref. |9|, which suggests similarities between 
phenomenological theories for magnons and electrons. 

We now relate the currents j and } q as well as the 
time derivative of the collective spin-density direction, m, 
to the thermodynamic conjugates via kinetic coefficients. 
The kinetic coefficients can be further identified by noting 
that the currents j and } q are determined by the chem- 
ical potential and temperature gradients as well as the 
magnetic texture dynamics, which exerts fictitious Berry 
phase gauge fields on the magnons. By assuming the spin- 
rotational symmetry of the magnetic texture and isotropy 
in real space, we obtain the magnon and energy current 
gradient expansion: 

-dip =T ik j k + U ik d k T/T -p(mx dim + /3<%m) • m , 

=Uj k j k - n ik d k T - pi (m x ^m + ftc^m) • m , 

5(1 + amx)m + m x H e ff = p [dim + /3m x dim] ji 

+ Pl [dim + ftm x dim] d{I '/T , (4) 

where the equation for the magnon current has been in- 
verted, Tik, and n ik are the resistivity, Peltier and 
thermal conductivity tensors, respectively, which are in 
general temperature and texture dependent {pifii has to 
be treated as one coefficient when p\ = 0) and the LLG 
part has been completed by employing the Onsager reci- 
procity principle. The heat current in Eq. (|3| includes 
the contribution from the magnetic energy contained in 
the texture j* (e.g., exchange energy flow accompanying 
the DW motion). In situations when the texture energy 
contribution can be separated from the "real" heat current, 
we choose to subtract this contribution, i.e., j£ = } q —} q in 
Eq. Q. By invoking the time reversal argument, one can 
show that the gradient expansion of the texture energy 
contribution has the form (j q )i = p t fi t dim • m. Equation 
(p3| can now be reformulated using the "pure" heat current 
jqf= i q — iq where as it follows from Eq. (p| the effective 
field has to be redefined as H e ff = H — ptfrdimdi (1/T). 
The new notations allow for more natural separation of 
the magnetic and heat degrees of freedom, e.g., in the 
simplest approximation one can assume that even in an 
out-of-equilibrium situation, when diT ^ and dip =^ 0, 
H e ff depends only on the instantaneous texture m(r, t). In 
general, we may expand H e ff phenomenologically in terms 
of small diT and dip. The form of Eqs. Q is the same for 
both sets of variables (jj H j g , H e ff <o> H) but the kinetic 
coefficients are different. The rest of the paper relies on 
equations written for (j£, H e ff). Note that in Ref. |9 it 
was implied that the latter set of variables had been used. 
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The tensors T^, and 

ac^/c can depend on tempera- 
ture and texture, i.e., to the leading order, as 

«ifc(n*fc,Tifc) =S ik T) + ^(n,T)(^ m ) 2 ] 

+ ^l(n,T)^ m • d k*n + &«(n,T)m • (c^m x d k m) . 

Equations Q should also be applicable to magnetization 
dynamics in the presence of charge currents (the magnon 
current has to be replaced by the charge current) as dis- 
cussed in Ref. |9|. These equations should describe such 
fictitious-field induced effects on magnons (electrons) as 
the Hall effect (coefficient by), the Ettingshausen effect 
(coefficient frn/ft)> the Nernst effect (coefficient bu/T), and 
the Righi-Leduc effect (coefficient &n/n). 

The gradient expansion in Eq. Q assumes short 
magnon wavelength compared to the characteristic tex- 
tures length scale. In YIG, for example, the former is 
~ 1 nm at room temperature, which should not pose a 
serious constraint for such adiabatic description. 

Bottom up construction of phenomenology for 
thermal magnons. — Since we concentrate on the low 
temperature limit (on the scale set by the Curie tempera- 
ture) the LLG phenomenology remains reliable microscop- 
ically. Coarse graining thermal fluctuations thereof would 
however lead to modified effective quantities (such as spin 
density and stiffness constant). Consider a ferromagnet 
with space- and time-dependent spin density sra(r, t) with 
magnitude s saturated at a constant value and direction 
described by a unit vector m(r,t). The spin density is 
assumed to have two components - fast and slow where 
the slow component slowly varies in space and time with 
much larger characteristic scales compared to the fast com- 
ponent. The effect of topological gauge fields due to mag- 
netic textures of the slow component can be captured by 
considering the Lagrangian [20] : 

C = j d 3 r [D(m) • m - £(m, d a m)} , (5) 

where D(m) = s[n x m]/(l + m-n) is the vector potential 
of the Wess-Zumino action with an arbitrary n pointing 
along the Dirac string, E(m, d a m) is the magnetic energy 
density describing the exchange energy as well as the ex- 
ternal and anisotropy fields. Equation ([5| leads to the 
Landau-Lifshitz (LL) equation sm — m x 5 m E — 0. We 
assume in our description that the magnetic energy den- 
sity is given by E/M s = A(d a m) 2 — m • H m /2 — m H 
where A is the exchange stiffness, M s is the saturation 
magnetization, H m describes magnetostatic and magne- 
tocrystalline anisotropics (largely ignored in our descrip- 
tion of thermal magnons at sufficiently high temperatures; 
but would otherwise suppress the fictitious forces acting 
on magnons [21]), and H is the external magnetic field. 
For the purpose of deriving the equations of motion, we 
use a coordinate transformation after which the z— axis 
points along the spin density of the slow dynamics. In the 
new coordinate system, small excitations will only have 
m x and m y components. 



Equations for the slow dynamics in the original (lab) 
frame can be obtained from the LL equation with a dis- 
sipative term by coarse-graining over fast variables, i.e., 
s (m) fast = s(r, t)m s with the final result: 

$m s + m s x H| ff =hd a }a + «ni s x (m s x H| ff ) 

- h(a - /3)m s x (daia) , (6) 

where the l.h.s. is written for the magnon- averaged 
spin density s, m and m s are unit vectors, i.e., m/ = 
m — ($/s)m s corresponds to the fast dynamics, the 
spin current density hj a = M s A[m.f x S a mj], H| ff = 
— 5 ms E{(s/s)m s , d a [(s / s)m s ]} is the effective field of the 
slow dynamics and j3 accounts for the fact that magnons 
misalign with the direction of the slow dynamics El. Note 
that Eq. ([6| assumes short magnon wavelengths compared 
to the slow texture length scale. By multiplying Eq. ([6| 
by 1 + am s x from the left and using an approximation 
daia ~ (jidi) m s (ji corresponds to the magnon number 
flux) which is true for slowly varying textures, we recover 
the Landau-Lifshitz-Gilbert (LLG) equation: 

5(l+am s x)m s +m s xHg ff = ft[l + /3m s x] (j i S i )m s , (7) 

In order to describe small excitations (spin waves) 
in the coordinates with the z— axis pointing along the 
slow dynamics, we introduce 3x3 rotation matrix R = 
ex.-p(i/jJ z )ex.-p(QJ y )ex.-p((j)J z ) with J a being the 3x3 ma- 
trix describing infinitesimal rotation along the axis with 
index a. In the new coordinates, we have m — >• m = Rm 
and <9 M —> (<9 M — A^) with A^ = (d^lVjRr 1 (the index 
ji = 0, .., 3 denotes the time and space coordinates). Since 
the matrix A^ is skew-symmetric, we can introduce a vec- 
tor v4 M so that A/^m = A,^ x m. In a specific gauge 
with the Euler angle ip = 0, the elements of v4 M become 
v4 M = (— sin 0(9^0, (9^0, cos Od^cj)). The equation describ- 
ing spin waves follows from the LL equation subject to the 
coordinate transformation: 

i(d t - ag)m+ = A (dji - A z a f m + + V(v)m + . (8) 

Here V(r) = m s • H/s - A(A 2 X + A 2 y )/2 is the effec- 
tive potential and m s is the unit vector along the spin 
density of the slow dynamics. Second order terms in 
the effective potential ~ A 2 ^ are not treated systemat- 
ically as there are similar corrections that lead to the 
coupling between the circular components of spin wave 
m± = m x (r,i) ± im y (r,t) where m^^(r,t) describes the 
transverse excitations in the transformed coordinates with 
the z axis pointing along the direction of m s (in the 
absence of texture m + ~ exp(iqr + iuj q t)). We omit- 
ted damping terms for the case when the region of in- 
terest is smaller than the length corresponding to the 
lifetime of spin waves. The coupling between the cir- 
cular components of spin wave due to anisotropics is 
disregarded since we assume that the exchange effects 
dominate. By quantizing spin waves, we introduce the 
field operator ip = ^s/2h 2 m+ = \A/^Xlqfrq e2qr corre- 
sponding to ra + and describing magnons with spectrum 
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u q = V(r)+Aq 2 where and 6 q are creation and annihila- 
tion operators. In such notations, the current is written as 
j a = —ih 2 A(^d a i/j—iljd a ^). Note that Eq. ^ describes 
charged particles moving in the fictitious electric S a = 
-d t A%-d a (A% + V(r)) = hm s -(d t m s xd a m s )-hd a V(r) 
and magnetic Bi = (h/2)e^ k m s • (c?fcm s x djih s ) fields pro- 
duced by the magnetic texture. In the diffusive regime, the 
transport of magnons due to such fields can be found by 
solving the Boltzmann equation within the relaxation-time 
approximation. 

The texture independent transport coefficients can be 
expressed through the following integrals: 



7' 



oc/3 



deT(e)(e - /i) r 



dfo 

' de 



! 



dS, 



v a vp 



(9) 

where \i is the chemical potential of magnons, r(e) is the 
relaxation time, e(q) = huj q , v a = duo q /dq ai dS £ is the 
area d 2 q corresponding to the constant energy e(q) = e 
and fo = {exp [(s — ^/fc^T] — l} -1 is the Bose-Einstein 
equilibrium distribution. Thus the resistivity tensor is 
p = J§ , the analogue of the Peltier coefficient for 
magnons is II = —JiJq 1 , and the thermal conductivity 
is k = (J2 — JiJq 1 Ji) /T . The relaxation time r con- 
tains the conserving number of magnons part r c , e.g., cor- 
responding to magnon-magnon interactions and disorder 
scattering, and the nonconserving part, e.g., correspond- 
ing to the Gilbert damping r a . We can sum up different 
contributions with different dependence on temperature 
according to 1/r = l/r c + l/r a . In this work, we assume 
that the scattering is dominated by the Gilbert damping 
contribution r a ~ (2acj) _1 , where a is the Gilbert damp- 
ing, which is consistent with our bottom up approach (in 
some cases, however, the contribution corresponding to 
the temperature independent scattering length gives bet- 
ter agreement to the experiment |22 1). By taking the 
quadratic spectrum for YIG at room temperature and 
a ~ 10 -3 1 23 1, we obtain the mean free path ~ 100 nm 
from the r a relaxation time for magnons. 

The above description of magnons contains dissipative 
/3 corrections |4 related to a small magnon mistracking 
of the slow magnetic texture. Such corrections have been 
extensively studied for charge flows in transition metals 
and one can use the analogy between the electronic and 
magnonic systems in order to introduce such corrections 
into the equations for currents. By following this prescrip- 
tion, one can arrive at Eqs. Q and ^ with the transport 
coefficients given by the Boltzmann equation. The coef- 
ficient P can then be estimated as f3 ~ a based on the 
absence of any additional energy scales, apart from the 
thermal energy, in the low-temperature (compared to the 
Curie temperature) limit captured by our bottom up ap- 
proach. More careful estimates of j3 can be obtained mi- 
croscopically by solving equation for nonequilibrium spin 
polarization in the presence of the transverse relaxation 
[6] or by using the scattering matrix approach [24] . 

By relating the microscopic description based on Eq. 



m(x,t) 



/ -> > ^ 

Tl 



Fig. 2: (Color online) A domain- wall moves towards the hotter 
end of the wire due to the flow of magnons j. Here, we consider 
transverse head-to-head Neel domain wall parallel to the y axis 
in the easy xy plane. The constants K and K± describe the 
easy axis and easy plane anisotropy. 

([7]) with Eqs. Q, we can immediately identify some of the 
kinetic coefficients in Eq. Q, i.e., pi = 0, P1P1 = 0, p = 
—h and a/P ~ 1. Note that p\ and p\P\ can be nonzero in 
systems with more than one magnon band and/or added 
anisotropics. In order to find the remaining coefficients, 
one could use Boltzmann equation and the relaxation time 
approximation applied to three dimensional magnons with 
Bose distribution [e.g., Eq. ([9])]. Such description can be 
further improved by considering the dissipative corrections 
described by 77 terms. 

Magnonic domain wall motion by temperature 
gradient. — In this section we show how Eqs. Q can be 
used to describe the domain wall motion by employing the 
Walker ansatz. To this end, we will describe the domain 
wall in Fig. [2] by the Walker ansatz valid for weak field 
and current /temperature biases [lQ|[25]: 



ip(r,t) = 3>(t), In tan 



0(r,t) 



X(t) 



W{t) 



(10) 



where the position-dependent spherical angles ip and 
parametrize the magnetic configuration as m(x) = 
m(x)(cos#, sin ^ cos (/p, sin#sin(/?), X(t) parametrizes the 
net displacement of the wall along the x axis, and we 
assume that the driving forces (H, j and j q ) are not 
too strong so that the wall preserves its shape and 
only its width W(t) and out-of-plane tilt angle un- 
dergo small changes. By substituting the ansatz (10) 
in Eq. ^ with the effective field given by H e ff = 
75 [(H + (5/ 's)Km x )x — (s/s)K±m z z + (5/s)iV 2 m] , we 
obtain: 



X 

w 



aX 
W 

- a<& 
W 



r y(s/s)K± sin 2$ pj 
2 sW ' 



A 



K + K_l sin 2 $ 



(11) 



where 7 is the gyromagnetic ratio, A is the stiffness con- 
stant, and K and Kj_ describe the easy axis and easy plane 



p-4 



Thermomagnonic spin transfer and Peltier effects in insulating magnets 



anisotropies, respectively. Note that the factor s/s before 
the stiffness constant appears due to coarse-graining over 
fast variables in Eq. (|6| . This dependence of the effective 
stiffness constant on the average spin density can lead to 
additional contribution to the domain wall velocity 1 30 . In 
this work, we concentrate on the low temperature limit (on 
the scale set by the Curie temperature) in which case s ~ s 
and the domain wall motion is dominated by the spin- 
transfer torque contribution. The steady state solution of 
Eq. (11) below the Walker breakdown with = const 
and X = vt leads to the result v = (jHW + h/3j/s)/a or 
equivalent ly 



7 W 



H 



F f3 
67r 2 Asa 



d x T, 



(12) 



where A is the thermal magnon wavelength and we 
introduce a numerical dimensionless factor Fo(x) = 
f dee 3 / 2 e e + x / (e e + x — l) 2 ~ 1 evaluated at the magnon gap, 
x = huo/kBT, which corresponds to dimesionless part of 
the integral J\ in Eq. From Eq. (12), we estimate 



the domain wall velocity in YIG at room temperature to 
be 1 cm/sec for the temperature gradient 1 K/fim. 

Peltier effect in quasi-one-dimensional wire. — 

Peltier effect describes the heat flow accompanying the 
flow of carriers. Since the moving magnetic texture will 
induce the flow of magnons, one can discuss the heat cur- 
rents resulting from such a process. In this section, we 
concentrate on quasi-one-dimensional systems with a DW 
propagating in a one dimensional wire connected to two 
reservoirs in quasi-equilibrium state (Fig. [2]). Contrary to 
the previous sections, the wire can be in the ballistic as 
well as in the diffusive regime. By writing the equation 
for the entropy production in the form analogous to the 
microscopic form in Eq. (J3|: 



L 
T 



TL 



Ti ./i 2 -/ii -2MH 

J — = X 



L 



L 



(13) 



and identifying the thermodynamics variables jg : j and 
X, we can phenomenologically generalize Eq. (U| to sys- 
tems possessing domain wall (solitonic) solutions describ- 
able by a single generalized coordinate X. The general 
phenomenological equations for the domain wall dynam- 
ics become: 

X = -O x (2M s H/L) - Xj £ - XT j q , 
St IT = G T j q - Tj £ - Oxt(2M s H/L) , 

j = Oj£ + Tj j q + Xj (2M s H/L) , (14) 

where L is the length of the wire, £ = — (^2 — Vi)/L, 
£t = —(T2 — Ti)/L and the kinetic coefficients now cor- 
respond to the whole wire, i.e., Ot = + k,ZT), Oj = 
1/(T + TZT) and G Tj = 1/(11 + II/ZT) with the conven- 
tional figure of merit relation for magnons ZT = U 2 /TTk. 
Using Eq. (14) one can fully describe the interplay be- 



in the quasi-one-dimensional systems. The kinetic coeffi- 
cients can be extracted from the bulk values in the diffu- 
sive regime (see below) or can be calculated by scattering 
theory methods in the ballistic regime |24|. 

In the following, we analyze the efficiency of heat pump- 
ing by a moving domain wall in a system depicted in Fig. 
[2] Given that the rate of dissipation due to the DW motion 
in Fig. [5] is 2M S HX and it can be assumed to be divided 
equally between reservoirs, we can calculate the ratio be- 
tween the useful heat taken from the cooled reservoir and 
the dissipated heat as the domain wall moves from the 
left end of the wire to the right end. Such ratios are often 
calculated in order to characterize thermoelectric circuits 
[26] . By maximizing the rate of cooling as a function of 
the domain wall velocity we obtain: 



COP cool 



~cold 
Jq 



T c VI + TZ mc - T h /T c 



2HMX T h - T c VI + TZ n 



(15) 

Here j^ old is the heat current leaving the cooled reservoir 
and we define the magnetocaloritronic figure of merit |9 
by analogy to the thermoelectric figure of merit ZT: 



TZ n 



O 2 

T O x ' 



(16) 



where we assume that £ = and such definition of TZ mc 
ensures that Eq. (15) is identical to the expression for 
the thermolelectric COP COO \ after TZ mc is replaced by 
ZT. The figure of merit in Eq. (16) is also related to 



the maximum efficiency of the magnetocaloritronic power 
generator in a device driven by ac magnetic field |9 con- 
trary to the geometry optimized magnetothermopower 
in applied dc magnetic field |27 . By taking some par- 
ticular DW solution we can relate Eqs. Q with Eqs. 
( |l4| ), e.g., for a transverse head-to- head Neel DW solu- 
tion in Eq. jiob, O x = LW/2as + T(l + ZT)Q 2 Xj , 
Xj (l + ZT) = (p/3/T - p&ZT/Ity/as and O xt (1 + 
ZT) = (pPZT/U + p 1 (3 1 /TK)/as. 

We now estimate TZ mc with p±/3i =0 according to Eq. 
([7]) and assuming that the scattering is dominated by non- 
conserving mechanisms described by r a ~ (2au)~ 1 . The 
corresponding mean-free path is assumed to be smaller 
than the system size in order to ensure the diffusive limit 
[in principle the ballistic regime can also be treated by 



Eq. ( 14 ) in which the kinetic coefficients should be found 



using an appropriate microscopic approach]. By invoking 
Eq. ([9| we express the figure of merit as follows: 



TZ n 



2F 1 VkFT(P/a) 2 (p/a) 2 



37r 2 yfA/hWLs W 2 Xs/h 



(17) 



tween the domain wall motion and magnon/heat currents 



where A is the thermal magnon wavelength, the 
wire length is taken to be L ~ W, and we in- 
troduce a numerical dimensionless factor F\{x) = 
[Ji/(Jo J2)} fdee 3 / 2 e e + x /[(e e + x - l) 2 (e + x)] evaluated at 
the magnon gap, x = Huo/kBT, which corresponds to 
dimesionless part of the ratio J 2 / J2 of integrals in Eq. 
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([9]). We also assume that the domain wall size can be 
estimated as W ~ ^ A/ujq where ujq can be, e.g., the de- 
magnetizing energy. By taking material parameters for 
YIG at room temperature |28 we arrive at W < 100 nm 
and TZ mc ~ 10 _4 (/3/a) 2 , which is quite low. However, 
larger ratios /3/a could be expected as one approaches 
the Curie temperature. Materials with smaller DW size 
should be more efficient in heat pumping according to 
Eq. (17), e.g., we estimate that in EuO W ~ 1 nm and 
TZ mc ~ 10" 3 (/3/a) 2 at - 10 K using the following pa- 
rameters: the localized spin S = 7/2, the lattice spacing 
ao = 5.1 A and the exchange integral Jq/Ub ~ IK |29 
(Jo = aosA/AS 2 for a face-centered lattice). In traditional 
thermoelectrics ZT plummets to zero much faster than 
TZ mc oc \fT in Eq. (17), making magnonic heat pumps 



promising for cryogenic applications. Furthermore, TZ mc 
scales as s" 1 with the spin density thus the dilute mag- 
netic systems (with sufficiently narrow W) should also be 
suitable for such applications. 

To conclude, we developed a phenomenological theory 
describing magnon and heat currents and the magneti- 
zation texture dynamics. Under some simple model as- 
sumptions, we are able to extract information about all 
the phenomenological parameters from the Gilbert damp- 
ing, the exchange integral, the localized spin and the lat- 
tice spacing. The j3 viscous coupling also appears in our 
description and is related to the magnon dephasing time. 
Our estimates show that the viscous coupling effects be- 
tween magnetization dynamics and magnon flows can be 
strong in materials with low spin densities and narrow do- 
main walls, which can allow the magnonic manipulation 
of magnetization dynamics and heat pumping. This opens 
new prospects for thermomagnonic devices, e.g., thermo- 
magnonic heat pumps and generators, that at low temper- 
atures could effectively compete with traditional thermo- 
electrics. When resubmitting our manuscript, we became 
aware of two recent works [30| |3l] that numerically study 
the magnonic spin-transfer torque and domain wall motion 
and arrive at results that are consistent with our studies. 
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